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INTRODUCTION 


Pressure fluctuations of turbulent flow are known to be the source of driving force for structural 
vibrations and radiation of acoustic energy. Theoretically, exact solution of the pressure fluctuations can 
be obtained by directly integrating the Poisson’s equation for pressure. The spectral properties of the 
pressure fluctuations can then be computed by Fourier transforming the two point correlation of the 
pressure fluctuations. The transformed equation involves terms with up to fourth-order velocity correla- 
tions. Approximate assumptions are employed to solve the higher-order correlation terms. 

Application of this approach to a constant-mean-shear homogeneous turbulent flow is described 
in detail by William K. George, et al. [1]. For wall pressure fluctuations on the surface of a flat plate, 
R. H. Kraichnan [2] estimated that the ratio of the root-mean-square (rms) wall pressure fluctuations to 
the wall shear stress is 6. Experimental measurements, reviewed by W. W. Willmarth [3], and computa- 
tions done by Panton and Linebarger [4] show that this ratio is lower than Kraichnan’s estimation. 
Experimental data of wall pressure fluctuations including compressibility effect have been collected and 
correlated by Langanelli, Martellucci, and Shaw [5]. They conclude that a value between 1.7 and 3 for 
the ratio of the rms pressure fluctuations to the wall shear stress may be appropriate for the subsonic 
turbulent boundary layer. Results of a direct numerical simulation of R. A. Handler, et al. [6] for a 2-D 
channel flow show a value around 3.32 for this ratio. 

The main drawback of the approach using the direct solution to the Poisson’s equation for 
pressure is that the geometry of the boundary has to be simple so that the Green’s function for integra- 
tion can be handled easily. Generally, the Green’s function for integrating the Poisson’s equation must 
satisfy the following two conditions [7,8] : 


V^G = 5(Ar) in the flow field . 


— = 0 on the solid boundary . 


where 5(Ar) denotes the Dirac delta function based on a displacement in space At, and n is in the direc- 
tion normal to the solid boundary. To find such a function for a complicated boundary geometry, such 
as that of the interior of the piping installations of the space shuttle main engine, would be practically 
impossible. 

In this paper, however, a completely different approach is taken. This new approach stems from 
the momentum equation for incompressible flow. The Reynolds decomposition concept is employed to 
derive a transport equation for the root-mean-square pressure fluctuations. This transport equation is then 
closed by approximately modeling the higher order terms which results in three empirical constants to be 
determined from experimental measurements. 


A sensitivity study on the empirical constants is included for the case of an incompressible high 
Reynolds number 2-D channel flow from which one of the predictions is compared with experimental 
measurements and other predictions. Results of an annulus flow prediction are also included. 

Major advantages of the present approach are that the concern on the geometric complexity is 
eliminated and the model can be solved using any CFD code along with the equations governing the 
flow field. 


DERIVATION OF THE TRANSPORT EQUATION FOR THE ROOT-MEAN-SQUARE 

PRESSURE FLUCTUATIONS 


Basic equations involved in deriving the transport equation for the root-mean-square fluctuations 
are the continuity and momentum equations of incompressible flow. They are 


Continuity: 


3U| 

3x ; 


= 0 


Momentum: 


Wi I ap , , a2tJ i 

dt j 3xj p dxj 3 xj3xj 
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where Uj represents velocity components, p is the fluid density, and P denotes the fluid pressure force. 

Using Renolds’ decomposition technique, the velocity components and pressure can be decom- 
posed into a mean quantity and a time-dependent component. They are written in the following form: 


U i = 


Ui + Ui 


P = p + p 


Where the variables with over bar denote time-averaged mean values. Substituting these relations into 
equation (2), a time-dependent momentum equation can be obtained. The result is 


9uj 3uj duj 3uj 3uj { 9p 3 u; 

"at” + Uj 3x^ + U J 3xj~ + dxj" ” U j 3xj p 3)q + V 3^ 


(3) 


Multiplying equation (3) by p, taking time-averaging, and taking the sum for i - 1,2,3, the follow- 
ing equation is obtained: 
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Using the time-dependent continuity equation, the following transport equation for the rms pressure 
fluctuations is obtained: 



The terms in equation (5) have been labeled with Roman numerals which permit us to state the relation 
in words: The rate of change of pressure fluctuations through pressure strain redistribution (I) is equal 
to (II) the convective diffusion and viscous diffusion of mean-square pressure fluctuations, plus (III) 
the production of pressure fluctuation, plus (IV) the turbulent viscous dissipation. 

Closure modeling is required to solve this equation. This is described in the next section where 
three empirical constants are introduced. 


MODELING THE TRANSPORT EQUATION 


The terms in equation (5) cannot be solved analytically. Approximate closure modeling to these 
terms is necessary to make equation (5) solvable. 

First, the terms on the left hand side of equation (5) can be written as: 


9uj 9uj 9puj 

^ 9t ^ ^ 9xj 9t 


9p u i 


9p 


+ U ; — U: - U: lb 

J 9xj 1 9t J 1 


dp 
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Assuming that the first two terms on the right hand side of equation (6) are small and negligible, equa- 
tion (6) is approximately modeled as 


9uj 9uj / 9 ■J'p 2 

P — — + u: p - — ss C i \/K I — 

9t •> 9xj 1 \ 9t 



(7) 


where Cj is called the coherence function in statistics and K is the turbulent kinetic energy. 

The diffusion terms (II) in equation (5) are modeled by dimensional analysis which takes the 
following form: 
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( 8 ) 



where C 2 is an empirical constant used to adjust the coefficient of diffusivity and £ represents a length 
scale which is related to the turbulent kinetic energy, K, and its dissipation rate, e, through the following 
equation: 


£ = C 


3/4 K 3 / 2 




where the empirical constant - 0.09. 

The production term (III) in equation (5) is expressed in terms of the mean quantities of turbu- 
lent flow using dimensional analysis. This gives 


(HI) «C 3 M tv ^ 


( 5 )' • 


(9) 


where C 3 is an empirical constant and n t represents the turbulent viscosity. p t can be calculated from 
the following equation: 


K" 


M t = P S — 


The dissipation term (IV) in equation (5) is also approximated using dimensional analysis. The 


result is 


(IV) * - c 4 




( 10 ) 


where C 4 is another empirical constant. 

Finally, by letting P' = the model of equation (5) can be written as: 


3P' _ 3P' _ 1 9 /£P' 9P'\ / 9 u i \ 2 
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where a p = C 4 C 2 / 2 , = C 3 /C 4 and C a = C 4 /CJ are the three empirical constants of the present model. 

Boundary conditions for solving equation (11) are: 

a) Production term » dissipation term near the wall. 

b) Gradient of P' across the line of symmetry is zero. 

Two of the empirical constants, C a and C^, can be estimated from data of homogeneous turbu- 
lence and wall pressure fluctuations measurements. 

For nearly isotropic homogeneous turbulence, equation (11) can be written as: 



Integrating the above equation once, 



Also, the K-e model [9] can be written as: 
dK 

dT = - e (13a) 


de _ e 3 

di " " Ce2 T 


(13b) 


Integrating equation (13b) once, the following relation is obtained. 



dt = - 


In (e/e Q ) 


C e 2 


C 3 / 4 


(14) 


From numerical and experimental investigations [10,11], the decay of isotropic turbulent kinetic energy 
can be represented by: 
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k = K o (t-t 0 r n 


or 


§■ - - nK 0 (t-t 0 )-<" + ') 


= - e 


(15) 


where the decay constant n « 1 .087 for very high Reynolds number flow. 

By letting e Q = nK Q in equation (15), the following expression can be obtained: 


In 


(c)" 


(n+1) ln(t-t 0 ) 


(16) 


Substituting equations (14) and (16) into equation (12), the result is: 


c a = - 


c e2 C„ 3/4 ta (p’o) 


(n+1) ln(t-t 0 ) 


(17) 


It is. also known [ 1 0, 1 1 ] that the pressure fluctuations in isotropic turbulence has approximately 
the same decay constant as the turbulent kinetic energy. Therefore, 


h (£)*-" 


ln(t-t Q ) 


and 


C 


a 


C,2 C u 3 / 4 n . 3/4 

(n+1) " 


0.1643 


where C e 2 = 1.92 in standard K-e model. 

The constant can also be estimated by applying the wall boundary condition to the present 
model. That is, 


C b ~ C a 


Vk p' 


L ^t® 




near the wall 
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where du/dy = C p ^\/lc/jcy by using the logarithmic law of wall with the von Karman constant k = 0.4. 

It has been concluded from experimental measurements of wall pressure fluctuations in boundary 
layer flows [3,5] that the range of rms pressure fluctuations, normalized by wall shear stress, is between 
1.7 and 3.0. 

This is equivalent to conclude that the ratio of wall pressure fluctuations to pK just outside the 
viscous sublayer is roughly between 0.5 and 0.85. If the higher value of this ratio is accepted, then 


near the wall » 0.85 


Therefore, only one empirical constant, a p , remained to be determined by matching the predic- 
tion to other sources of data in the flow field. To the author’s knowledge, no pressure fluctuations 
measurements have been made away from boundaries in any confined flow region. However, results of 
higher level numerical simulations, e.g., large eddy simulations for channel and annulus flows done by 
Schumann [12], can be used to estimate the constant o p . 



SENSITIVITY STUDY ON a p AND DATA COMPARISONS 


Effects of <7p on te predictions of rms pressure fluctuations are investigated numerically for an 
incompressible fully developed 2-D channel flow. Equations for solving this problem are: 


Momentum equation: 


Id/.- du \ _ dP _ 2 

y j dy y Mt dy / dx pUr 

(18a) 

K - e model: 


1 d / dK\ _ / du\ 2 

y i dy \ a K dy / * Mt \dy) 

(18b) 

/£\ /duV 

y i dy y o € dy/ e \K / M \dy/ 

(18c) 

Pressure fluctuations model: 
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( 18 d) 



where = 1.0, a £ = 1.3, C e 2 = 1.92, C e j = 1.43, C a = 0.1643 and = 0.85. The parameter j = 0 is 
for 2-D channel flow problems and j = 1 is for axisymmetric pipe flow problems. 

A finite difference routine is used to solve the above equations. Equations (18a), (18b) and 
(18c) are solved first until a converged solution is obtained. Equation (18d) is then used to solve for the 
rms pressure fluctuations. 

Experimental measured [13] and numerically computed results of the mean velocity profile, 
turbulent kinetic energy and its dissipation rate for the fully developed 2-D channel flow are illustrated in 
Figure 1. Reynolds number based on the maximum mean velocity and the channel half width is 62,000 
which is about the same as the experiment. 

Effects of Op on the predictions of rms pressure fluctuations are shown in Figure 2. Notice that 
the variation of ffp does not affect the predictions near the wall where zero diffusion conditions are 

applied. The near-wall rms pressure fluctuations are fixed by the ratio of C a and C^. The levels of 

predicted p' increase slightly in the vicinity of the wall with increasing value of a p while the levels of p' 
decrease strongly with increasing <7p elsewhere. 

One of the cases studied above is selected (o p = 6.0) and compared with experimental measure- 
ments and other predictions. This is shown in Figure 3 where the wall pressure fluctuations measured by 
Ball [14] and Elliott [15] and predicted by Handler [6] are used for comparisons. Pressure fluctuations 
data comparisons away from the wall are provided by Schumann [12] and estimated by a constant- 
mean-shear isotropic turbulence formula given by George [ 1 ] . This formula is 

P' = p 0.33 ^ P. 2 u 2 + 0.42 u 4 (19) 

where u 2 = 1/3 (u 2 + v 2 + w 2 ) = 2/3 K, £ = u 2 /e and du/dy stands for the mean shear (i.e., an average 
slope of the mean velocity profile). In the present data comparisons, the local mean velocity gradient is 
used for du/dy to estimate p' across the channel. 

It is shown in Figure 3 that the predictions of the present model with Op = 6.0 compare well 

with experimental measurements and other predictions. Especially, the shape of the p' profile is very 
similar to the predictions of Schumann using large eddy simulation for high Reynolds number flow. This 
suggests that Op = 6.0 is a good approximation for the present model. 

This set of empirical constants (i.e., Op = 6.0, C a = 0.1643 and = 0.85) of the present model 

are used to predict pressure fluctuations in a fully-developed annulus flow and compared with 
Schumann’s prediction. This is shown in Figure 4. Ratio of the larger pipe radius to the smaller pipe 
radius of the annulus is 5.0. It is shown in Figure 4 that agreement between the present predictions and 
the Schumann’s predictions is generally good except the discrepancies near the wall. 
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SUMMARY 


The transport equation for the root-mean-square pressure fluctuations has been derived from the 
time-dependent momentum equation. Approximate closure modeling for this transport equation has 
been proposed with three empirical constants. 

Two of the empirical constants have been estimated from data of homogeneous turbulence and 
wall pressure fluctuations measurements. The third constant for the diffusion coefficient of the present 
model has been estimated by comparing the results of a large eddy simulation for 2-D channel flow. 

Data comparisons between the current predictions and other sources of data have been included 
for a 2-D channel flow and an annulus flow. Results of the comparisons show that the predictions of 
the present model are at least qualitatively good. More tests of the present model are required to fine- 
tune the constants of the present model before it can be confidently used in engineering design analysis. 
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Figure 1. Velocity profile, turbulent kinetic energy and its dissipation 
rate of a fully-developed 2-D channel flow. 


Ca = 0.1643 


Cb = 0.85 


Figure 2. Effect of Op on the predictions of rms pressure fluctuations 
for the 2-D channel flow. 
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